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Abstract 

We investigate on the relative inclination of the planets B and C 
orbiting the pulsar PSR B1257+12. First, we show that the third Ke- 
pler law does represent an adequate model for the orbital periods P 
of the planets, because other Newtonian and Einsteinian corrections 
are orders of magnitude smaller than the accuracy in measuring Pb/c- 
Then, on the basis of available timing data, we determine the ratio 
sintc/sinzB = 0.92 ± 0.05 of the orbital inclinations ?b and ic in- 
dependently of the pulsar's mass M . It turns out that coplanarity 
of the orbits of B and C would imply a violation of the equivalence 
principle. Adopting a pulsar mass range 1 < M < 3, in solar masses 
(supported by present-day theoretical and observational bounds for 
pulsar's masses), both face-on and edge-on orbital configurations for 
the orbits of the two planets are ruled out; the acceptable inclinations 
for B span the range 36 deg < ib J$ 66 deg, with a corresponding 
relative inclination range 6 deg < (ic — «b) ^ 13 deg. 
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The 6.2-ms PSR B1257+12 pulsar was discovered in 1990 during a high 
Galactic latitude search for millisecond pulsars with the Arecibo radiotele- 
scope at 430 Hz [14] . Two years later, PSR B1257+12 turned out to be 
orbited by at least two Earth-sized planets-B and C-along almost circular 
paths [15]. In 1994 it was announced the discovery of a third, Moon-sized 
planet- A-in an inner, circular orbit [IB]. Its presence, questioned by Scherer 
et al. [IT] , was subsequently confirmed in [6j[I7]. The relevant orbital pa- 
rameters of the PSR B1257+12 system are listed in Table [TJ Note that the 
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Table 1: Relevant parameters [8] of the three planets |X5|, I16| A, B and C, 
hosted in the PSR B1257+12 system p3], derived from analysis of timing 
data ranging 12 years (1990-2003) collected at the 305-m Arecibo telescope. 
P is the orbital period, x is the projected barycentric semimajor axis of 
pulsar's motion and 7 = m/M [7] is the ratio of the planet's mass to the 
pulsar's mass. Figures in parentheses are the formal la uncertainties in the 
last digits quoted. 



Planet 


P(d) 


x (ms) 


7 (io-°) 


A 


25.262(3) 


0.0030(1) 




B 


66.5419(1) 


1.3106(1) 


9.2(4) 


C 


98.2114(2) 


1.4134(2) 


8.3(4) 



ratios 7b, c of the masses of B and C to M were measured from timing data 
exploiting their mutual gravitational perturbations [7], without using the 
standard reference value M = M = 1.4M© for the pulsar's mass. 

Given the peculiarity of the PSR B1257+12 system, it is certainly im- 
portant to deepen the knowledge of the orbital configuration of its planets 
in order to gain insights about the evolutionary dynamics of such a rare sys- 
tem. Here, without making any a priori assumptions about the inclinations 
of B and C, we wish to constrain them from the available timing data by 
assuming a reasonable interval of masses for the pulsar. 

To this aim, we will exploit the third Kepler law whose use is justified in 
detail below. First of all, let us note that, by defining s = s'mi, it is possible 
to write the planetary relative (i.e. pulsar-to-planet) semimajor axis a in 
terms of the measured quantities 7 and x as 

i+ MW = /l + 7W ) 
m ) s \ 7 / s 

where c is the speed of light. Note that both x and 7 were phenomenolog- 
ically determined in [8] independently of the third Kepler law itself and of 
the pulsar's mass. The third Kepler's law is 



P\ 2 G 3 



- 



2tt J GM(1 + 7) ' 



(2) 



from which it is possible to express the pulsar's mass in terms of the phe- 
nomenologically determined quantities x, 7, P, apart from s which will be 
considered as unknown. 
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Note that a purely Keplerian model for the orbital period is quite ade- 
quate because non-Keplerian corrections like those due to the oblateness of 
the pulsar (if any), to the planet-planet interaction and to the 1PN 0{c~ 2 ) 
terms are negligible given the present-day accuracy in determining P% /q . In- 
deed, concerning the oblateness of the central mass, its contribution AP( obl ) 
to the orbital period of an orbiting test particle can be written as [5] 

Ap( obl) = _6tt£J2 
VGMa 

where \J-z\ < 1 represents the first even zonal harmonic coefficient of the 
multipolar expansion of the gravitational potential of the pulsai0 and R 
is the pulsar's radius; by assuming typical values R = R = 10 km and 
M = M = 1.4M© one has 

AP^ « - J 2 x 10- 12 d. (4) 

The corrections AP 1 - 3 ^ bod y) to the orbital period of a planet of mass m 
induced by another planet of mass m can be written as [5] 

A p(3rd body) _ 47rGm ^ ( 5 ) 

na 

where n = Gm/a? is the Keplerian mean motion of the perturbed planet 
and a is the perturber's semimajor axis. In the case of the planets B and 
C it turns out that 

Ap (3rd body) _ 1() -15 d ^ 

The fact that the planets B and C are in a 3:2 resonance does not affect their 
orbital periods. Indeed, according to the Lagrange's perturbation equation 
for the variation of a [T], 

aoc^-, (7) 

where a = —nT p is related to the time of pericentre's passage T p and H\ 
is the interacting Hamiltonian of eq.(23) in ref. [7]. Since H\ does not 
explicitly contain a there is no secular change in the semimajor axis of B 
and C. The mutual perturbing effects employed in ref. [7] to estimate 7b 
and 7c are not the corrections to the Keplerian orbital periods. 



x Of course, in this rough order-of- magnitude estimate there is no need to take into 
account contributions due to the peculiarity of the matter state in the neutron star. 
2 As a consequence, we also used is = 53 deg and ic = 47 deg obtained in [8] with M. 
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The 1PN Post-Newtonian correction AP( 1PN ) to the orbital period of 
order 0(c~ 2 ) can be written as [2] 

A pdPN) = ^l^GM~ a . (8) 



it turns out that 



AP£™ ] « !0 -6 d. (9) 

The latter contribution is the most important post-Keplerian correction to 
P, but it is two orders of magnitude smaller than the la formal errors in 
the phenomenologically determined orbital periods quoted in Table [TJ 
Thus, we have 

XC 



GM 



2vr(l+ 7 ) n2 



(10) 



Taking the ratio of eq. (|10p for both B and C allows us to get information 
about the relative orbit inclination independently of M itself: indeed, we 
have, from Table [1] 

S - - = ( £ ( ^) ( 1 - ± ^) ^ = 0-92 ± 0.05. (11) 



«B \PcJ \XBlcJ V 1 + 7B/ 

The la error was conservatively assessed by propagating through eq. (fTTj) 
the uncertainties in Pq , Pq , xb , xq , 7b > 7c quoted in Table Q] and linearly 
adding the resulting biased terms. It turns out that the most important 
sources of errors are 7b and 7c yielding 5S 7B =0.02 and 5S 1C = 0.03. 

It maybe interesting to note that, by assuming siniB = sinic, the quan- 
tity S in eq. (Illh may be interpreted as a measure of a violation of the 
equivalence principle. Indeed, according to [ID], by putting m g = m\{l + rj) 
for the gravitational and inertial masses of a test particle orbiting a cen- 
tral body of mass M, the 2— body problem encompassing a violation of the 
equivalence principle can be described by the same formulas derived in the 
classical 2— body problem, provided that whenever they contain the product 
GM, we substitute it with GM(1 + rf). Applying it to eq. (fTPj) written for 
the planets B and C, it turns out that S can be interpreted as 

V 1 + VB ) 

Thus, coplanarity and the result of eq. (jlip would yield a violation of the 
equivalence principle in the PSR B1257+12 system at 1.6a level. However, 
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it must be noted that such a test would be stronger if one had evidence that 
the composition of the two planets was very different. 

The constrain of eq. (|lip is a consistent, genuine dynamical one which 
does not make use of any assumption about M; the authors of [8] did not 
obtain it. However, we note that the values of 7b/c> entering eq. (fTTj) . were 
measured in [8] by using the model of [7] which neglects terms in sin 2 (//2), 
where I = |«b — *c| is the relative inclination of the orbital planes of B and 
C assumed to be I < 10 deg. In Figure [1] it is plotted the allowed region 
in the plane {iB,ic}, according to eq. (fTTj) . delimited by the minimum and 
maximum values of the ratio S. The blue dashed line represents the copla- 
narity condition, while the green line is the maximum value of I allowed by 
the system's parameters. It can be noted that for face-on (i — > deg) ge- 
ometries the orbital planes tend to be coplanar. Most remarkable deviations 
from coplanarity (more than 10 deg with a maximum of about 32 deg for 
ib = 90 deg) occur for edge-on (i — > 90 deg) geometries, but in such cases 
caution is required since we would fall outside the / < 10 deg condition on 
which our analysis relies upon; Figure [1] shows that this occurs for is > 56 
deg. In the following we will explore the viability of various possible orbital 
inclinations in terms of physically plausible values of M. 

We will study the behavior of eq. (fTUf) for B and C as a function of is 
(because of eq. (jlip ) in order to constrain the inclinations. Indeed, since 
we have no independent information at all on ^b/C! we wm use a reason- 
able interval of masses for the pulsar to constrain them (and their relative 
inclination / through eq. (jlip which is independent of M). We wish to 
preliminarily notice that, in principle, PSR B1257+12, as a member of the 
rare class of the planetary pulsars, may have had a different formation and 
evolution with respect to the other neutron stars. However, in absence 
of any other indication on the details for the evolutionary history of PSR 
B1257+12, we will rely the following analysis upon standard mass intervals, 
commonly adopted for other kinds of neutron stars. This is a notable differ- 
ence with respect to [8] in which the pulsar's mass was kept fixed to 1.4Mq. 
Theoretically speaking, different Equations-Of-State for nuclear matter in- 
side neutron star yield different pulsar's mass ranges; we will adopt 1 — 3Mq 
[9]. Let us note that present-day observations are all compatible with that 
range. As to the lower bound, all the best determinations of the mass of a 
neutron star fall well above 1 M , approaching that value only for the most 
uncertain cases (see ref. [12] for an overview of measured pulsar masses). 
As to the upper bound, the highest securely measured value of the mass 
of a pulsar is that recently obtained for PSR J04374715 (about 1.76 ± 0.20 
M [E]). Anyway, our adopted mass range also includes the cases of more 
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Figure 1: Plot of ic = arcsin(S'siniB) for the minimum (lower red curve) 
and maximum (upper red curve) values of S according to eq. (fTTl) . Such 
constraints are independent of the pulsar's mass. The blue dashed line rep- 
resents the coplanarity case. The distance between the blue coplanarity line 
and the lower red curve yields the maximum value of the relative inclination 
/ allowed by the system's parameters along with their uncertainties; it is 
depicted in green and shows that I > 10 deg for is > 56 deg. 
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Figure 2: Allowed region for i-Q as a function of M. The 1 — 3M interval 
of [9] has been chosen. 

massive neutron stars (e.g. ref. [3]), whose claimed high mass must still be 
confirmed by additional investigations. 

The lex error in the value of M calculated using eq. (jlOl) can be con- 
servatively evaluated by propagating the uncertainties in P,x,j of Table 
[H and linearly summing the resulting individual biased terms. It is of the 
order of 13% and the major contribution to it turns out to be due to 7. By 
considering the allowed regions for M determined by the curves M ± 5M by 
means of eq. (fTUI) applied to both B and C, it turns out that the tightest 
constraints come from B whose allowed region is entirely enclosed in that 
due to C. In Figure [2] we depict the constraints on i-q for 1 < M < 3 in 
solar masses. As can be noted, it turns out that 36 < i-Q < 66 deg. From an 
inspection of Figured] it turns out that the relative inclination / is different 
from zero being 6 < I < 13 deg. It is interesting to note that larger values 
of /, which, at least to a certain extent, may still be compatible with the 
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analysis presented herd_|, are ruled out by the lower bound on the pulsar's 
mass. 

The authors of [8], by using M = M, obtain 49 deg < i-Q < 57 deg and 
44 deg < iq < 50 deg. On one hand, our analysis confirms-as expected-that 
a larger range of values for zb and ic is allowed when a suitable interval of 
values for the mass of the pulsar is taken into account. On another hand, 
thanks to the new constrain implied by eq. f 1 1 [) . it is now possible to show 
that not all the combinations of orbital inclinations i-Q and ic indicated in 
ref. [8] are acceptable (even when adopting the same uncertainty level of 
ref. [8]). In fact, as depicted in Figure O the rectangle of acceptable values 
(49 deg < ib < 57 deg) x (44 deg < ic < 50 deg) indicated in ref. [8] is not 
entirely included in the region allowed by the mass-independent constrain 
ofeq. (pi]). 
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